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Abstract 

The  theory  of  high-gain  error-actuated  feedback  con¬ 
trol,  developed  by  Porter  and  Bradshaw,  was  applied  to  the 
design  of  various  longitudinal  decoupling  flight  control 
systems  for  an  advanced  aircraft.  This  aircraft  incor¬ 
porates  both  pitching  moment  (elevators)  and  direct  lift 
(flaps)  control  devices,  permitting  the  application  of 
multivariable  control  theory.  The  controllers  developed 
in  this  study  utilized  output  feedback  with  proportional 
plus  integral  control  to  produce  desirable  closed-loop 
response  with  minimal  interaction  between  outputs. 

Because  of  the  structure  of  the  system,  measurement 
variables  different  from  the  outputs  are  necessary  to 
apply  this  method.  -This  report  describes  how  entire 
eigenstructure  assignment  can  be  used  to  determine  appro¬ 
priate  measurement  equations  by  assigning  their  correspond¬ 
ing  transmission  zeros.  A  singular  value  decomposition 
was  used  to  choose  the  eigenvectors  from  their  permissible 
subspaces.  The  results  show  that  these  modes  appear  in 
the  output  response,  even  for  very  high  gain.  Therefore, 
selection  of  good  eigenvalues/eigenvectors  was  shown  to 
be  crucial  to  the  successful  application  of  this  theory. 

*  "*  '  At- 

'  In  addition,  this  report / examined  the  effect  of  varying 
the  other  design  parameters  on  the  closed-loop  system 
response.  . 
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USE  OF  ENTIRE  EIGENSTRUCTURE  ASSIGNMENT  WITH 
HIGH-GAIN  ERROR-ACTUATED  FLIGHT  CONTROL  SYSTEMS 

I  Introduction 

In  the  past  several  years,  a  great  deal  of  work  has 
been  done  with  Control  Configured  Vehicles  (CCV) .  A  CCV 
is  one  in  which  advanced  control  technology  as  well  as 
aerodynamics,  structures,  and  propulsion  are  employed  in 
the  initial  definition  process  (Ref  1:1).  Advanced  con¬ 
trol  technology,  such  as  direct  force  control,  can  permit 
new  and  highly  advantageous  ways  in  which  an  aircraft 
maneuvers.  Specifically,  it  is  the  direct  force  control 
capability  of  these  aircraft  that  is  addressed  in  this 
report.  Direct  force  control,  in  conjunction  with  conven¬ 
tional  elevator  control  surfaces,  allows  independent  con¬ 
trol  of  the  normal  force  and  pitching  moment  degrees  of 
freedom  of  the  aircraft.  The  most  recent  version  of  this 
type  of  aircraft  is  the  AFTI-16  (Advanced  Fighter  Tech¬ 
nology  Integration) ,  which  is  an  F-16  aircraft  modified 
to  include  additional  control  surfaces,  such  as  canards, 
leading-,  and  trailing-edge  flaps. 

Background 

The  addition  of  these  control  surfaces  allows  a 
large  number  of  direct  force  control  flight  modes,  each 
of  which  involves  different  coupling  between  the  aircraft 
degrees  of  freedom.  For  example,  in  the  longitudinal  axis, 
the  modes  found  to  be  most  useful  are  pitch  pointing, 
vertical  translation,  and  direct  lift  control.  The  pitch 
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pointing  mode  is  effected  by  varying  the  aircraft  pitch 
attitude  while  holding  the  flight  path  angle  constant. 

The  vertical  translation  maneuver  requires  varying  the 
flight  path  angle  while  holding  the  pitch  attitude  con¬ 
stant.  Varying  flight  path  angle  while  holding  the  angle 
of  attack  constant  causes  the  aircraft  to  exhibit  the 
direct  lift  control  maneuver.  These  flight  modes  are 
quite  different  than  conventional  aircraft  modes  in  that 
they  require  non-interaction  between  specific  outputs. 

To  accomplish  these  unusual  maneuvers,  a  controller  must 
be  designed  which  causes  one  output  variable  to  track  a 
given  command  while  eliminating  the  response  in  other  out¬ 
put  variables.  This  controller,  of  necessity,  must  gen¬ 
erate  signals  to  deflect  the  two  longitudinal  control 
surfaces  (flaps  and  elevator)  in  order  to  produce  the 
desired  output  response.  That  is,  this  controller  must 
be  capable  of  controlling  multiple-input,  multiple-output 
(MIMO)  systems. 

Numerous  design  methods  are  currently  in  use  for 
flight  controllers.  Classical  control  techniques  are  not 
well  suited  to  this  type  of  design,  as  they  are  very 
difficult  to  use  with  MIMO  systems.  Two  major  modern 
control  design  methods  are  optimal  control  and  entire 
eigenstructure  assignment.  In  the  optimal  control  design 
method  (Ref  12)  ,  it  is  desired  to  minimize  the  control 
effort  required  to  track  given  commands.  This  requires 
minimizing  the  performance  index 


J  =  ^-[y(tf)  -  r(tf)]T  H[y(tf)  -  r(tf)] 

+  \  /Qtf  Uy(t)  -  r(t)]T  Q(t)(y(t)  -  r (t)  ] 

+  uT(t)R(t)u(t) }dt  (1-1) 

for  the  system  described  by 

x  =  Ax  +  Bu 

y  =  Cx  (1-2) 

where  y(t)  is  the  output  vector  and  r(t)  is  the  desired 
value  of  the  output  vector.  H,  Q(t),  and  R(t)  are  sym¬ 
metric,  positive  definite  quadratic  weighting  matrices. 

The  optimal  control  law  which  minimizes  J  is 

u*(t)  =  K ( t) x ( t)  -  R-1 ( t) BT  d (t)  (1-3) 

where  the  time-varying  feedback  gain  matrix  is 

K (t)  =  -R-1(t)BT  P(t)  (1-4) 

and  P(t)  is  the  solution  of  a  non-linear,  time-varying 
matrix  differential  equation  known  as  a  Eiccati  equation, 
defined  by 

P(t)  =  -ATP(t)  -  P(t)A  +  P ( t ) BR  _1 ( t ) BTP { t )  -  CTQ(t)C 

P(tf)  =  CTHC  (1-5) 

Because  the  boundary  condition  in  equation  (1-5)  is  speci¬ 
fied  at  tj,  the  Riccati  equation  must  be  integrated  back¬ 
ward  in  time,  thus  making  on-line  computer  implementation 
difficult.  In  equation  (1-3),  d(t)  is  found  from 
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(1-6) 


-d(t)  =  ATd ( t )  -  K(t)BR_1  (t)BTd(t)  -  CTQ(t)C 
d(tf)  =  -CTHCr(tf) 

Once  again,  this  equation  must  be  integrated  backward  in 
time  due  to  the  given  boundary  condition,  thus  requiring 
information  to  be  known  about  the  commands  r(t)  in  the 
future.  In  addition,  since  t^  should  not  be  much  larger 
than  the  initial  time  (usually  on  the  order  of  five  seconds 
for  aircraft  maneuvers) ,  a  steady-state  analysis  is  not 
theoretically  in  order,  and  the  time-varying  Riccati 
equation  must  be  integrated.  Solutions  to  this  type  of 
equation  are  extremely  costly  computationally  and  often 
difficult  to  achieve.  Lastly,  selection  of  the  three 
weighting  matrices  must  be  made,  which  is  the  "art"  of 
the  method.  All  of  these  factors  contribute  to  the  dif¬ 
ficulty  and  complexity  of  an  optimal  tracking  control 
problem. 

Entire  eigenstructure  assignment  is  based  upon 
selection  of  the  system  closed  loop  eigenvalues  and 
corresponding  eigenvectors.  Given  a  system  described  by 
the  same  state  and  output  equations  as  before  (equation 
(1-2))  a  control  law 


u  =  Kx  +  u  (1-7) 

—  —  — c 

is  assumed.  Substituting  this  control  law  into  the  state 
equations  yields  the  closed  loop  system 

x  =  ( A+BK) X  +  Bu^  (1-8) 


4 


Using  the  gain  matrix 


K  =  [  w  ^  * 1  1  I  £^  r  £2  9  *  *  *  9  S-n  ^  ( i  —  9 ) 

produces  a  closed  loop  system  (equation  (1-8) )  with  eigen¬ 
values  X^  and  corresponding  eigenvectors  5^,  i  =  l,2,...,n, 
provided  the  eigenvectors  and  vectors  ok  satisfy  the 
relationship  (Ref  6) 


e  X/{  [A-XiI,B]  } 


(1-10) 


Care  must  be  taken  in  selecting  the  eigenvectors  from  the 
null  space  [A-X^IjB]  so  that  the  matrix  [£2 ,£2 »•••»£  1  is 
invertible.  Furthermore,  these  eigenvalues  and  eigenvec¬ 
tors  must  be  chosen  so  that  the  amount  of  interaction 
between  outputs  is  minimized  through  adjusting  the  tran¬ 
sient  responses.  Some  work  has  been  done  (Ref  2)  on 
developing  methods  to  select  the  eigenvectors  to  promote 
desirable  response,  but  the  procedure  at  this  point  in 
time  is  still  largely  trial  and  error. 

In  this  report,  a  method  designed  especially  to 
track  commands  and  requiring  relatively  simple  calculations 
(easily  performed  on-line)  is  examined.  This  method  should 
greatly  aid  in  the  development  of  flight  controllers  for 
CCV  type  aircraft. 

Problem  Statement  and  Approach 

It  is  desired  to  design  and  evaluate  a  flight  con¬ 
troller  (using  a  method  proposed  by  Porter  and  Bradshaw 
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(Ref  3) )  for  an  AFTI-16  aircraft  which  will  exhibit  the 


maneuvers  of  pitch  pointing,  vertical  translation,  and 
direct  lift  control.  The  relatively  new  method  of  con¬ 
troller  design  developed  by  Porter  utilizes  singular  per¬ 
turbation  methods  in  the  design  of  systems  incorporating 
high-gain  error-actuated  controllers.  The  goal  of  this 
study  is  to  minimize  the  amount  of  "tuning"  required  to 
effect  the  three  direct  force  control  maneuvers  and  to 
examine  closely  the  results  of  varying  parameters  in 
Porter's  method.  To  minimize  the  "art"  involved  in  apply¬ 
ing  the  design,  the  technique  of  entire  eigenstructure 
assignment  was  used  to  select  the  closed-loop  transmission 
zeros.  The  permissible  subspace  for  selecting  the  eigen¬ 
vectors  was  obtained  by  using  a  singular  value  decomposi¬ 
tion  (Ref  6) . 

The  successful  achievement  of  these  maneuvers  is 
based  on  several  criteria.  First,  the  closed-loop  system 
must  be  stable.  Secondly,  all  transient  responses  should 
die  out  after  no  more  than  three  seconds  (that  is,  achiev¬ 
ing  98%  of  their  steady  state  values) .  Another  require¬ 
ment  is  that  overshoots  of  the  steady  state  values  are 
minimal,  and  the  closed-loop,  steady-state  responses  should 
achieve  their  commanded  values  with  no  observable  error. 
Lastly,  the  closed-loop  system  should  be  very  robust, 
thereby  conforming  to  the  above  requirements  when  variations 
are  made  in  the  aircraft  plant. 


Scope  and  Assumptions 

The  purpose  of  this  study  was  to  thoroughly  evaluate 
the  design  method  proposed  by  Porter  for  CCV  applications. 
As  a  result,  design  iterations  were  primarily  to  examine 
the  effect  of  various  parameters  rather  than  to  iterate 
to  a  "best”  solution.  To  this  end,  only  the  longitudinal 
axis  and  its  corresponding  CCV  modes  were  examined.  The 
assumptions  listed  below  were  made  to  reduce  the  complexity 
of  the  model.  None  of  these  should  seriously  affect  the 
conclusions  reached  as  a  result  of  this  effort. 

1.  The  earth  is  flat  and  non-rotating 

2.  The  aircraft  is  rigid  and  of  constant  mass 

3.  The  atmosphere  is  at  rest  with  respect  to  the 
earth 

4.  Gravity  is  a  constant  acceleration 

5.  The  aircraft  X-Z  plane  is  a  plane  of  symmetry 

6.  The  motion  of  the  aircraft  about  a  straight 
and  level  trim  condition  can  be  adequately 
modeled  by  a  linear  dynamical  system  of  equations 

Organization 

First,  the  theory  used  in  the  design  is  described 
and  the  method  for  developing  a  controller  discussed. 

Next,  Chapter  3  develops  the  linearized  models  for  the 
aircraft.  Chapter  4  then  presents  the  results  of  applying 
this  method  to  the  aircraft  models.  Several  examples 
are  examined  and  the  results  described  in  detail, 
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especially  showing  the  effects  of  varying  design  parameters 
and  the  aircraft  plant.  Lastly,  Chapter  5  discusses  the 
results  and  provides  recommendations  for  further  studies. 
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II  Theoretical  Development 


Porter's  Method 

In  Porter's  method  (Ref  3),  high-gain  error-actuated 
analog  control  is  used  in  conjunction  with  singular  per¬ 
turbation  methods.  First,  attention  is  given  to  the 
rationale  behind  this  approach.  High  gain  frequently  has 
the  advantage  of  making  the  closed-loop  system  very  robust, 
or  insensitive  to  parameter  variations  within  the  plant. 
Also,  high  gain  may  be  desirable  from  the  viewpoint  of 
disturbance  rejection,  since  it  can  produce  a  very  tight 
tracking  loop,  although  this  disturbance  rejection  aspect 
was  not  examined  in  this  study. 

Error-actuated  control  is  the  keystone  of  this 
method.  The  control  law  developed  by  this  method  auto¬ 
matically  minimizes  the  difference  between  the  commands 
and  the  outputs  as  the  gain  is  increased.  Therefore, 
tracking  is  accomplished  when  the  gain  is  made  large 
enough,  and  should  become  increasingly  tight  and  non¬ 
interacting.  How  large  the  gain  must  be  for  good  perfor¬ 
mance  is  one  of  the  principal  questions  adddressed  by  this 
study . 

Singular  perturbation  methods  (Ref  4)  are  used  to 
exhibit  the  asymptotic  structure  of  the  transfer  function 
matrices.  The  method  of  singular  perturbations  can  be 
used  to  examine  the  overall  effect  of  letting  a  certain 
parameter  become  smaller  and  smaller.  By  regarding  the 
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controller  forward  path  gain  (g)  as  the  reciprocal  of  the 
perturbation  parameter,  c,  the  results  of  singular  pertur¬ 
bation  theory  can  be  applied,  since  as  c  -*  0,  g  =  1/e  -*■  «. 
Using  this  theory,  the  asymptotic  properties  of  the 
closed-loop  transfer  functions  can  be  examined  to  deter¬ 
mine  the  location  of  the  closed-loop  poles  as  the  gain 
gets  very  large.  These  asymptotic  closed-loop  poles  then 
determine  the  system  response. 

The  basic  plant  for  Porter's  method  may  be  repre¬ 
sented  by  state  and  output  equations  of  the  respective 
forms 


and 


x2  (t) 


A12 

A22.  .-2(t) 


Z(t)  =  [c1  c2] 


*1  (t) 

x2  (t) 


(2-1) 


(2-2) 


where  B2  is  a  square,  non-singular  matrix.  Any  system,  by 
proper  state  and/or  control  transormations  can  be  put  into 
this  form.  For  the  examples  of  this  study,  such  transfor¬ 
mations  were  not  needed  since  the  original  system  of 
equations  was  already  in  this  form.  This  method,  as 
originally  proposed  by  Porter  and  Bradshaw  (Ref  5) , 
requires  that  the  matrix  C2B2  have  full  rank.  A  later 
extension,  which  is  examined  in  this  study,  compensates 
for  a  rank  deficient  C2B2  by  introducing  new  measurement 
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equations.  In  all  of  the  cases  examined  in  this  report, 
the  rank  of  C2B2  is  not  full,  so  that  extra  plant  measure¬ 
ments  of  the  form 


w(t)  =  [C1+MA11  C  2+MA1 2 ] x ( t ) 


(2-3) 


-  [F1  F2> 


(fc) 


*-x2  <  t)-1 


(2-4) 


were  necessary.  The  matrix  M  is  chosen  by  the  designer. 
The  high-gain  error-actuated  controller  is  governed  by 
a  control-law  equation  of  the  form 


u(t)  =  g{K  e(t)  +  K.z(t)}  (2-5) 

—  o—  1— 

where  K  and  Kn  are  controller  gain  matrices,  e(t)  is 
o  1  — 

the  error  vector  between  input  commands  v(t)  and  the 
measurements  w(t) ,  and  z(t)  is  the  time  integral  of  e(t) . 

It  is  desired  that  the  input  vector  u(t)  cause  the 
output  vector  y.(t)  to  track  any  constant  command  input 
vector  v(t) ,  so  that 

lim  { v  ( t )  -  y;<t)  }  =  0  (2-6) 

t-+oo 

Such  tracking  is  possible  using  the  above  controller  form 
(equation  (2-5)),  which  includes  integral  control ,  due  to 
the  fact  that  the  error  vector  e(t)  assumes  the  steady- 
state  value  (for  arbitrary  initial  conditions) 
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lim  e(t)  =  lim  {v(t)  -  w(t) } 

£-»-CO  t"^°° 

=  lim  (v(t)  -  ^(t)  -  M(A.  ,  x.  (t)  +  A12x2(t))} 
t"*00 


=  lim  (v(t)  -  y_(t)  }  =  0 

t->00 

since  (A^x^(t)  +  A^2— 2^^  9oes  to  zero  in  steady  state  as 
a  result  of  equation  (2-1).  In  equations  (2-1),  (2-2),  (2-3), 

(2-4),  and  (2-5),  x-^tje  n"£,  x2(t)eR£,  u(t)eR£,  y(t)cR£, 
w(t)eR£,  A ll£R(n-£)x(n-£),  A12cR(n'£)x£,  A21cR£x(n“£), 
A22eR£x£,  B2eR£x£,  C-^R^ (n_£) 

pvf  0  0 

F2eR  ,  rank  C2B2£^'  rank  F2B2=£ ,  e  ( t)  eR  ,  v(t)cR  ,  z^(t)  = 
2(0)  +  /ote(T)dieR£,  Kq£R£x£,  Ki£:R£x£,  g>0,  and  MeR£x  (n-£  ]  . 


,  C2sR£x£,  FlCR£x(n-£), 


The  error  states  are  governed  by  the  differential 
equations 


z(t)  =  e(t)  =  v(t)  -  w(t)  =  v(t)  -  F1x1(t)  -  F2x2(t)  (2-8) 


Substituting  this  expression  for  e(t)  into  equation  (2-5) 
and  then  the  resulting  equation  into  (2-1),  along  with 
equation  (2-8) ,  yields  the  closed-loop  equations 


-ut)  * 

- 1 

o 

1 

1 

NJ 

_ 1 

-  z  ( t )  - 

*1  (t) 

= 

0  An  A12 

xx  (t) 

-gB2K1  A21_gB2KoFl  A22-gB2KoF2- 

.-2 (tl 

LgB2KoJ 


v  ( t ) 


(2-9) 
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and  output  equations 


y(t) 


[o  c1  c2) 


z(t) 
x1  (t) 

X 2  (t). 


(2-10) 


The  results  of  singular  perturbation  theory  as  des¬ 
cribed  by  Porter  and  Shenton  (Ref  4)  are  used  to  determine 
the  asymptotic  form  of  the  transfer  function  matrix,  G(X), 
as  g-+°°.  Letting  the  perturbation  parameter  c  =  1/g  go  to 
zero,  the  following  system  is  obtained 


’  B1  " 

A1 

A2 

A  = 

B  = 

_A3/e 

A4/e_ 

.^2/ e. 

C  =  [C^  c2l 


where  the  components  of  the  asymptotic  closed-loop  A 
matrix  are 


A1  = 


(2-11) 


‘°  -F,' 

A  = 

-F2 

„°  All. 

_A12_ 

-B2KoF1 

]  a4  = 

[-B2 

(2-12) 


the  B  matrix  components  are 


B1  = 


I£1 


LO 


B2  =  tB2Ko] 


(2-13) 


and  the  C  matrix  elements  are 


c1  =  [0  c x3  c2  =  [c2] 


(2-14) 
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t 

J' 

l 


As  the  transfer  function  matrix,  G(A) 

assumes  the  form 

r  ( x )  =  r  ( x ) 

where 

f  ( X  )  =  C  ( X I 
o  n 

f(X)  =  C2(AI£  - 

The  matrices  in  equation  (2-16)  are 


C[AI-A]-1B, 


+  f  (A) 

(2-15) 

-  A  )_1i 

(2-16) 

o  o 

-1 

<3a4  )  g  B2 

(2-17) 

A  =  A. -A_A ~ 1A,  = 
o  12  4  3 


-K_1K, 
o  1 


,A12F2  Ko  K1  A11~AJ2F2  FlJ 


:2-18) 


B  =  B1-A_A~1A,  = 
o  12  4  3 


A  F'1 
A1 2t  2 


and 


(2-19) 


CQ  =  61-62A41A3  =  [C2F-1K;1K1  c1-C2F2lpl]  (2-20) 

The  "slow"  modes  of  the  tracking  system  are  defined  as 
those  which  are  not  directly  dependent  on  g.  These  are 
the  poles  of  r ( X )  and  are  easily  found  from 


det  [ X I  -  A  ]  =  0 
n  o 


(2-21) 


since  Aq  is  block  lower  triangular.  Simplifying  this 
relation  yields  the  two  sets  of  slow  poles  described  by 


Z.  =  { AeC ;  | XK  +  K.  I  =  0} 
X  O  X 


(2-22) 
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and 


Z2  *  f XcC:  I  XIn-e'All  +  ''l2F2ll'li  *  01  (2'23) 

The  "fast"  diodes  of  the  system  (those  directly  dependent 
on  g)  are  found  from 

det  [XI  -  gA4]  =  0  (2-24) 

so  that 

Z3  =  {XcC:jXI£  +  gF2B2KQ!  =  0}  (2-25) 

In  order  for  tracking  to  occur,  the  closed- loop 
system  must  obviously  be  stable,  requiring  that  Z^,  Z2» 
and  Z3  have  all  poles  in  the  left-half  plane.  Therefore, 
for  sufficiently  large  gains,  the  controller  and  trans¬ 
ducer  matrices  Kq/  K^,  and  M  must  be  chosen  so  that  Z^eC, 

2 2eC  ,  and  Z^C  ,  and  the  chosen  M  matrix  must  simul¬ 
taneously  satisfy  the  condition 


rank  F2B2  =  ^ 


(2-26) 


As  there  is  no  real  physical  insight  to  the  best  selection 
of  M  or  even  values  of  M  that  will  satisfy  the  above  con¬ 
ditions,  an  extension  to  Porter's  method  was  developed. 


Entire  Eigenstructure  Assignment  Extension 

To  eliminate  the  necessity  of  the  designer  selecting 


the  transducer  matrix,  M,  while  simultaneously  satisfying 


equation  (2-26)  and  ensuring  that  the  poles  of  Z2  are  in 
the  left-half  plane,  an  extension  to  Porter's  method 
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equation  (2-30)  becomes 


tXiI_All  A12] 


ii 


to . 

—l 


=  0 


(2-32) 


Therefore,  lies  in  the  null  space  of  [X^I-A^  A-^l  * 

Using  the  method  of  singular  value  decomposition  suggested 
by  Silverthorn  and  Reid  (Ref  6)  to  find  the  null  space  of 
Ai2^'  ^i  an<^  — i  can  calculated  for  any  choice 
of  closed-loop  eigenvalues.  Finally,  the  gain  matrix  S 
can  be  found  from  repeated  application  of  equations  (2-30) , 
(2-31),  and  (2-32)  for  the  n  closed-loop  eigenvalues. 

That  is 

S  —  I ,  .  .  •  ,  ui ^ ]  I ^  ,  C_2  >  •  •  •  > J  (2—33) 

From  equation  (2-33),  it  is  clear  that  the  matrix 
.  .  .  ,£nl  must  be  invertible,  therefore  the  ’ s  i  =  1,2,..., 
n  must  be  chosen  so  that  they  are  independent.  Neverthe¬ 
less,  much  design  freedom  still  exists  in  choosing  these 
eigenvectors.  Some  effects  of  choosing  the  are  shown 
in  the  results  section. 

Having  found  S,  the  matrix  M  can  be  found  from 
equation  (2-28)  and  equation  (2-3).  By  assuming  that  F^^ 
is  nonsingular,  from  equation  (2-28) 
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(unfortunately,  the  method  of  entire  eigenstructure  assign¬ 
ment  does  not  assure  that  ( F 2 ^ ^  exists.)  Using  the 
expressions  for  F^  and  F 2  in  equation  (2-3) ,  this  becomes 

(C2  +  MA12)S  =  (Cx  +  MAU)  (2-35) 

This  is  rearranged  to  yield 

M  =  (C2S  -  C1)  (An  -  A12S)-1  (2-36) 

Since  the  latter  matrix  is  imply  A  T  given  by  equation 
(2-29) ,  it  is  guaranteed  to  be  invertible  as  long  as  no 
eigenvalues  are  chosen  to  be  zero  (choosing  zero  for  a 
closed-loop  pole  would  make  the  system  marginally  stable, 
so  this  is  not  a  problem) .  Now  and  F2  can  be  calculated 
from  equation  (2-3) .  F2  must  be  checked  to  insure  that 
it  is  nonsingular.  If  it  is  singular,  either  different 
eigenvectors  must  be  selected  or  this  method  is  not  appli¬ 
cable  for  the  given  system  (such  a  case  is  presented  in 
the  results  section) .  The  method  described  above  for 
computing  M  is  an  alternative  to  the  method  suggested  by 
Porter  (guessing  M  so  that  F2  has  full  rank)  since  there 
is  considerably  more  physical  insight  behind  the  selection 
of  these  closed-loop  eigenvalues  than  the  selection  of  a 
transducer  matrix  M. 

In  summary,  to  use  this  synthesis  method,  first 
closed-loop  "slow"  poles  are  chosen  and  their  correspond¬ 
ing  eigenvectors  are  calculated  through  a  computationally 
efficient  singular  value  decomposition  (equation  (2-32)). 


18 


A  gain  matrix  S  is  computed  (equation  (2-33)),  from  which 
M,  ,  and  F2  are  calculated  (equations  (2-36)  and  (2-3)). 

Next,  the  fast  modes  are  selected.  The  matrix 
F2B2Ko  ke  chosen  to  be  a  diagonal  matrix  in  order 

for  increasingly  non-interacting  tracking  to  occur  as 
g-*oo.  From  equation  (2-17)  ,  a  diagonal  F2B2Ko  obviously 
produces  a  diagonal  transfer  function  f(A).  Therefore, 
letting 

F2B2Ko  =  ^  =  diag  {o^  ,o2  ,  .  .  .  ,  }  (2-37) 

the  matrix  K  can  be  calculated  for  chosen  o.  >0,  i  =  1, 
o  1 

2,  ...  ,1  from 

Ko  =  (F2B2)_1  Z  (2-38) 

Finally,  the  first  set  of  "slow"  poles  can  be  easily 
chosen  by  letting  be  a  constant,  positive  multiple  of 
Kq,  so  that  the  "slow"  poles  become  the  negative  of  the 
multiplicative  factor.  This  is  shown  by  letting  =  bKQ 
in  equation  (2-21) ,  so  that 

| XI£  +  K"1 (bKQ)  |  =  | A I £  +  b|  =  0  (2-39) 

This  clearly  shows  that  A^  =  -b,  i  =  1,2 ,...,£. 

Now  that  M,  Kq,  and  have  been  calculated  and  all 
the  necessary  conditions  for  tracking  satisfied,  the  closed 
loop  matrices,  equations  (2-9)  and  (2-10)  can  be  calculated 
for  choices  of  the  gain  parameter,  g.  From  these  equations 


the  closed -loop  eigenvalues  for  the  entire  system  and  the 
time  histories  of  all  states  and  outputs  can  be  calculated. 
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Ill  Development  of  the  Aircraft  Equations  of  Motion 


Linearized  Equations 

As  stated  in  the  introduction,  only  the  longitudinal 
equations  of  motion  are  analyzed  here.  The  model  of  the 
airframe  is  a  conventional,  small-perturbation  set  of 
equations  linearized  about  straight  and  level  flight.  The 
stability  axis  is  used  in  this  study,  where  the  trim  flight 
path  angle  (pitch  angle  of  the  stability  axis),  G  =  0. 

The  resulting  equations,  in  the  LaPlace  domain,  are  shewn 
in  Fig.  3-1  (Ref  7:298).  These  equations  are  transformed 
to  the  time  domain  and  separated  so  that  they  are  of  the 
form 


x  =  Ax  +  Bu  (3-1) 

The  results  of  this  transformation  are  shown  in  Fig.  3-2. 

Values  for  all  the  terms  in  the  equations  are  pre¬ 
sented  in  Appendix  A.  The  basic,  nondimensional  deriva¬ 
tives  were  taken  from  Ref  8,  and  represented  a  flight 
condition  for  which  h  =  3,000  ft  and  M  =  0.60.  These 
values  are  substituted  into  the  equations  in  Fig.  3-2  and 
are  presented  in  Fig.  3-3.  The  open-loop  poles  for  this 
system  and  their  natural  frequencies  w  and  damping  ratio 
£  are  given  in  Table  3-1. 

Short  Period  Approximation 

For  several  cases  examined  in  the  next  chapter,  the 
short  period  approximation  is  used  to  show  basic  charac¬ 
teristics  of  the  controller  and  to  compare  results  with 
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Fig.  3-1.  Linearized,  Small-Perturbation  Equations  of 
Motion  in  Stability  Axis. 
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Fig.  3-3.  Linear  Model  for  Basic  Aircraft 


Fig.  3-4.  Linear  Model  for  Short  Period  Approximation 
with  a  Replacing  w 


Table  3-1.  Eigenvalues,  Natural  Frequencies,  and  Damping 
Ratios  for  System  Shown  in  Fig.  3-3. 


Eigenvalue 

Natural  Frequency 
(rad/sec) 

Damping  Ratio 

-0.00757  ±  0 . 068 j 

0.0684 

0.1106 

5.453 

— 

-- 

-7.663 

— 

those  from  an  example  in  a  paper  by  Porter  and  Bradshaw 
(Ref  9).  In  the  short  period  approximation,  the  forward 
perturbation  speed,  u,  is  assumed  to  be  zero.  This 
eliminates  the  phugoid  mode  and  replaces  it  with  one 
pole  at  the  origin.  Also,  the  variable  w  (the  vertical 
perturbation  speed)  is  transformed  to  a  (angle  of  attack) 
by  the  relationship 

a  =  w/U  (3-2) 

where  U  is  the  aircraft  velocity.  These  equations  are 
shown  in  Fig.  3-4  and  the  open- loop  poles,  natural  fre¬ 
quencies,  and  damping  ratios  as  shown  in  Table  3-2. 

Table  3-2.  Eigenvalues  for  the  System  Shown  in  Fig.  3-4 

Eigenvalue 

0 

5.453 


-7.663 


In  several  test  cases,  the  effects  of  control 
actuator  dynamics  were  evaluated.  By  modeling  these 
actuators,  choosing  values  similar  to  those  used  in 
fighter  aircraft,  the  controller  can  not  generate  instan¬ 
taneous  deflections  of  the  control  surfaces.  This  factor 
could  be  extremely  important  in  the  evaluation  of  high- 
gain  controllers.  The  simple  model  used  for  the  actuator 
dynamics  was 


6  _  20 

6  s+20 

c 


(3-3) 


In  the  time  domain,  for  the  elevator  and  flap  control  sur¬ 


faces,  the  linear  model  is 


(3-4) 


Letting  6  and  6  F  be  additional  aircraft  states,  and  6 
e  r 

and  be  the  system  inputs,  the  plant  state  equations 

*-c 

for  the  full  model  and  the  short  period  approximation  are 
shown  in  Figs.  3-5  and  3-6,  respectively.  The  eigenvalues 
for  these  plants  are  the  same  as  for  their  unaugmented 
counterparts,  except  two  poles  at  -20  are  added  in  each 
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Fig.  3-5.  Full  Model  Aircraft  State  Equations  with  Control 
Actuators 
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IV.  Discussion  of  Results 


System  A  -  Short  Period  Approximation  With  No  Actuators 
The  first  group  of  test  cases  presented  in  this 
study  deals  with  designing  controllers  for  the  short 
period  approximation  model  without  control  actuator 
dynamics  shown  in  Fig.  3-4.  It  is  easily  seen  that, 
for  these  cases. 
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U  0} 
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A22 
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-1.3411- 

[0  0] 
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In  order  to  effect  the  pitch  pointing  and  vertical  trans- 

T 

lation  modes,  the  output  matrix  for  the  output  =  [0  y]  is 


iCiiCjl 


1  J  0  0 

1  I  0  -1 


(4-3) 


To  effect  the  direct  lift  control  mode,  the  output  matrix 


[c1;c2j 


1  1  0  0 

i 

o ;  o  -i 


(4-4) 


T 

is  required  for  y  =  [y  a]  .  From  equations  (4-3)  and  (4-4) 
it  is  obvious  that  C2  is  rank  deficient  in  both  cases, 
so  that  extra  plant  measurements  are  necessary. 

Since  the  system  is  augmented  with  two  integrators 
in  order  to  achieve  the  desired  tracking  properties,  the 
composite  system  is  of  order  five.  Consequently,  there 
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are  five  closed-loop  eigenvalues.  In  all  of  these  test 
cases  only  one  slow  eigenvalue  must  be  chosen,  since  A^ 
is  simply  a  scalar.  Because  there  are  two  system  inputs, 
and  since  the  associated  closed-loop  eigenvector  is  a 
linear  combination  of  the  vectors  lying  in  the  null  space 
of  [A^I-A^,  A^2^  /  two  constant  multipliers  (called  eigen¬ 
vector  multipliers)  must  be  chosen  to  form  the  eigenvector. 
Each  test  case  in  this  group  is  evaluated  separately. 

Case  A-l 

In  this  case  (-4+jO)  is  selected  for  the  slow  eigen¬ 
value  for  the  sake  of  comparison  with  Porter's  example  in 
Ref  9.  The  system  output  equations  are  those  of  equation 
(4-3).  Fig.  4-1  shows  the  M,  F2,  F^,  Z,  Kq,  and 

matrices  (henceforth  called  the  system  matrices)  for  the 
choice  of  {1,0}  for  the  eigenvector  multipliers.  Also 
shown  in  Table  4-1  are  the  closed-loop  eigenvalue  locations 
as  the  gain  parameter,  g,  is  increased  (referred  to  as 
eigenvalue  movements) . 

Pitch  pointing  is  effected  with  the  command  vector 

v  =  [1  OJT  deg  (4-5) 

As  g  is  increased,  the  system  goes  from  instability  to 
stability  (as  one  of  the  open-loop  eigenvalues  is  in  the 
right-half  plane),  and  once  stable,  increasing  g  reduces 
the  settling  times  and  final  value  overshoots.  Further 
comment  on  the  effects  of  g  is  presented  in  the  conclusions 
section.  The  time  responses  for  g  =  100  are  shown  in  Fig. 
4-2.  From  this  figure,  it  is  seen  that  y^  = 
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1.  System  Matrices  for  Case  A-l,  Multipliers  {1,0} 
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Table  4-1.  Eigenvalue  Movements  for  Case  A-l,  Multipliers  {1,0} 
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its  commanded  value  quickly  with  little  overshoot,  but 
y 2  =  y  has  a  large  perturbation  and  is  not  held  to  the 
desired  zero  value.  Commanding 

v  =  [0  1] T  deg  (4-6) 

effects  the  vertical  translation  maneuver,  and  the  results 
for  g  =  100  are  shown  in  Fig.  4-3.  Here,  y1#  is  held  to 
its  commanded  value  of  zero  quite  well,  but  y 2  rises 
extremely  rapidly  and  then  undershoots  slightly  before 
achieving  steady  state. 

In  an  attempt  to  improve  the  response,  the  eigen¬ 
vector  was  changed.  As  this  system  is  quite  simple  it  is 
possible  to  calculate  the  required  eigenvector  multipliers 
to  achieve  the  M  matrix  used  by  Porter  (Ref  9) .  Using 
eigenvector  multipliers  {1,-0. 25}  yields  the  system 
matrices  and  eigenvalue  movements  shown  in  Fig.  4-4  and 

Table  4-2,  respectively.  The  resulting  time  histories  for 
T  T 

v  =  [1  0]  deg  and  v  =  {0  1}  deg  are  shown  in  Figs.  4-5 

and  4-6,  respectively,  for  g  =  100.  The  pitch  pointing 
maneuver  is  now  effected  quite  well,  but  the  vertical 
translation  mode  is  still  extremely  fast.  This  problem 
can  be  rectified  by  using  a  "ramped-up"  inL  .t  as  Porter 
did,  in  which  case  these  results  match  Porter's  quite  well. 

Case  A- 2 

Once  again  (-4+jO)  is  chosen  for  the  slow  eigenvalue 
and  in  this  case  the  output  equations  of  equation  (4-4) 
are  used  so  that  the  direct  lift  control  maneuver  is 
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Fig.  4-4.  System  Matrices  for  Case  A-l 
{1,-0. 25} 
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Fig.  4-5.  Time  Responses  for  Case 
{1,-0. 25} 
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Responses  for  Case  A-l,  Vertical  Translation  Maneuver  with  Multipliers 


rn 

simulated  by  using  the  command  vector  v  =  [1  0]  deg. 
Choosing  the  multipliers  {1,0}  yields  the  system  matrices 
shown  in  Fig.  4-7.  The  eigenvalue  movements  are  identical 
to  those  of  Table  4-1.  The  time  responses  are  shown  in 
Fig.  4-8  for  g  =  100,  and  clearly  show  the  maneuver  is 
satisfactorily  effected. 

Case  A-3 

In  this  case  (-4+jO)  is  the  chosen  eigenvalue  again, 
but  here  the  controller  matrix  is  chosen  to  be  twice 
the  Kq  matrix,  yielding  {  (-2+j0) ;  (-2+j0) }  for  the  asym¬ 

ptotic  values  for  the  first  set  of  slow  eigenvalues  rather 
than  { (-1+jO) ;  (-1+jO) }  as  in  the  two  previous  cases.  The 

output  equations  of  equation  (4-3)  are  used  with  the 
command  vectors  of  equations  (4-5)  and  (4-6)  .  Four  dif¬ 
ferent.  choices  of  eigenvectors  are  made  to  examine  the 
results  of  changing  eigenvectors.  The  choices  of  the 
multipliers  and  a  summary  of  their  corresponding  results 
are  shown  in  Table  4-3.  The  results  indicate  that  the 
larger  the  values  of  the  elements  in  the  matrix  F2B2' 
worse  the  response.  The  last  case  is  shown  in  detail  in 
Figs.  4-9  and  Table  4-4  and  its  corresponding  time 
responses  shown  in  Figs.  4-10  and  4-11  for  g  =  100. 

These  responses  are  very  similar  to  those  of  Case  A-l  for 
the  same  eigenvector  except  that  in  this  case  the  pitch 
angle  6  overshoots  slightly  less  and  the  settling  times 
are  slightly  smaller. 
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7.  System  Matrices  for  Case  A-2,  Multipliers  {1,0} 
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Table  4-3.  System  Matrices  and  Maximum  Interactions  for  Several  Choices  of  Eigenvectors 
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System  Matrices  for  Case  A-3,  Multipliers 
{1,-0. 25} 
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Responses  for  Case  A-3,  Pitch  Pointing  Maneuver  with  Multipliers 
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-11.  Time  Responses  for  Case  A-3,  Vertical  Translation  Maneuver  with  Multipliers 
{1,-0. 25} 


Case  A-4 


This  case  was  used  solely  for  comparison  with  Case 
A-2,  and  therefore  is  identical  except  for  the  ratio 
between  Kq  and  K^.  The  results  are  given  in  Fig.  4-12, 
Table  4-5,  and  Fig.  4-13.  There  is  no  notable  difference 
in  the  time  responses  of  this  case  versus  those  of  Case 
A-2. 


Case  A- 5 

In  this  test  case,  the  output  equations  given  by: 
[c1',c2]x(t)  =  [6  q]T  = 

are  used  to  show  an  interesting  aspect  of  this  method. 
Using  the  method  developed  in  this  report,  for  all  choices 
of  slow  eigenvalues  and  corresponding  eigenvectors  the 
rank  of  F2  is  never  full.  Expanding  the  theory  developed 
by  Porter  for  this  example  shows  the  reason  for  this. 

Using  the  matrices  given  in  equations  (4-1) ,  (4-2) ,  and 

(4-7)  with  the  M  matrix  represented  by 
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the  following  F^  and  F 2  matrices  are  found 
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4-12.  System  Matrices  for  Case  A-4,  Multipliers  {1,0} 
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Table  4-5.  Eigenvalue  Movements  for  Case  A-4,  Multipliers 


sponses  for  Case  A-4,  Direct  Lift  Control  Mane 
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which  clearly  shows  F 2  does  not  have  full  rank  for  any 
choice  of  M.  This  corresponds  to  the  fact  that  the 
system  matrix  (Ref  13) 


XI -A 
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oj 


(4-11) 


is  rank  deficient  for  \=0  (transmission  zero  at  the  origin). 
Therefore,  the  system  matrix  may  be  rewritten 


f  A  b" 

L~C  0. 


which,  for  a  choice  of  q  (pitch  rate)  as  an  output, 
does  not  have  full  rank.  Therefore,  the  system  is  not 
fully  controllable,  which  is  a  necessary  condition  for 
successful  application  of  this  theory.  The  same  result 
is  true  for  a  selection  of  normal  acceleration ,  A^,  as 
an  output. 


Case  A-6 

This  case  is  identical  with  Case  A-l  except  that  I 

is  chosen  to  be  diag  (0.1, 0.1}  rather  than  diag  {1.0, 1.0}. 

The  resulting  system  matrices  are  identical  to  Case  A-l 

except  K  and  K.  are  one-tenth  that  of  the  former  case, 
o  1 

If  the  gain  in  this  case  is  chosen  to  be  ten  times  that 
of  any  example  from  Case  A-l,  the  closed-loop  eigenvalues 
and  all  responses  are  identical,  indicating  no  beneficial 
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M 


effects  from  increasing  or  decreasing  E  by  a  constant  fac¬ 
tor.  This  result  was  found  to  be  true  in  all  cases. 


System  B  -  Short  Period  Approximation  VJith  Control 
Actuator  Dynamics 

The  next  set  of  test  cases  have  control  actuator 
dynamics  modeled  in  the  aircraft  state  equations.  The 
model  used  is  shown  in  Fig.  3-6,  from  which  the  matrices 
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are  easily  recognized.  The  output  matrix  for  pitch  point¬ 
ing  and  vertical  translation  is  given  by 


tc1;c2j  = 
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-1  i  0  0 


4-15) 


and  that  for  direct  lift  control  by 
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0  -110  0 


(4-16) 


In  both  cases  C2  is  rank  deficient  and  extra  plant  measure¬ 
ments  are  required. 
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In  these  test  cases,  three  slow  eigenvalues  must 
be  selected,  thereby  allowing  complex  as  well  as  real 
eigenvalues  to  be  chosen.  As  two  vectors  lie  in  the  null 
space  of  (X^I-A^,  A^l.  *  =  1 »  2 ,  3,  the  selection  of 
six  eigenvector  multipliers  is  required  to  form  the  three 
closed-loop  eigenvectors. 

Case  B-l 

In  this  case  the  chosen  eigenvalues  are  {(-4r3j); 
(-3+jO)},  and  the  output  equations  are  those  represented 
by  equation  (4-14).  Appendix  B  discusses  the  selection 
of  complex  eigenvalues.  In  the  cases  where  complex  eigen¬ 
values  are  selected,  the  last  two  of  the  four  corresponding 
vectors  lying  in  the  null  space  are  assigned  zero  multi¬ 
pliers  to  reduce  the  complexity  of  the  problem.  This 
approach  still  reproduces  the  requested  eigenvalues. 

Since  the  control  surface  deflections  are  now 
included  as  state  equations,  the  control  time  histories 
are  easy  to  evaluate.  In  the  previous  test  cases  without 
actuator  dynamics,  the  control  law  has  to  be  evaluated 
and  integrated  in  order  to  obtain  these  histories.  This 
is  but  one  advantage  to  modeling  control  actuator  dynamics. 

Several  choices  of  eigenvector  multipliers  are  made 
in  this  case.  First,  multipliers  (1,  10,  0,  0,  1,  1}  are 
chosen  and  the  resulting  system  matrices  and  eigenvalue 
movements  are  shown  in  Fig.  4-14  and  Table  4-6,  respectively. 
The  corresponding  control  histories  and  time  responses  for 
g  =  1000  are  shown  in  Figs.  4-15  through  4-18  for  command 
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Table  4-6.  Eigenvalue  Movements  for  Case  B-l,  Multipliers  {1,10 
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Fig.  4-18.  Time  Responses  for  Case  B-l,  Vertica] 
{1,10,0,0,1,1} 


r  — 

T  T 

vectors  v  =  [10]  deg  and  v  =  [01]  deg.  These  results 
show  the  control  surface  deflections  to  be  very  abrupt 
in  the  first  second  and  show  the  outputs  to  be  signifi¬ 
cantly  interactive  (especially  in  the  pitch  pointing  mode) . 

Next  the  multipliers  are  changed  to  {1,  1,  0,  0,  1, 

0}  with  corresponding  results  for  g  =  1000  shown  in  Fig. 
4-19,  Table  4-7,  and  Figs.  4-20  through  4-23.  These 
results  indicate  two  major  trends.  First,  the  magnitude 
of  the  elements  in  the  F^B  matrix  is  greatly  reduced 
from  that  of  the  previous  case.  Second,  the  control 
surface  deflections  are  not  as  abrupt  in  this  case.  The 
time  responses  in  this  case  are  greatly  improved  over 
the  previous  one. 

Changing  the  multipliers  slightly  to  {1,  1,  0,  0, 

1,  1}  produces  even  better  responses.  Results  in  this 
case  are  shown  in  Fig.  4-24,  Table  4-8,  and  Figs.  4-25 
through  4-28  for  g  =  1000.  The  magnitude  of  the  elements 
of  the  ^2^2  Inatr;'-X  show  little  change  over  the  previous 
case,  but  the  control  deflections  are  slightly  less 
abrupt.  Consequently,  the  time  response  is  made  much 
more  non-interactive. 

Case  B-2 

The  output  equations  given  by  equation  (4-15)  along 
with  the  choices  of  { (  —  4±3j) ;  (~3+j0)}  for  the  slow  eigen¬ 
values  are  used  in  this  case.  Numerous  choices  of  eigen¬ 
vector  multipliers  were  made,  with  two  examples  shown 
here.  In  the  first  example,  the  choice  of  { 1 , 0 , 0 , 0 , 1 , 0 } 
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4-19.  System  Matrices  for  Case  B-l,  Multipliers 
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Fig.  4-23.  Time  Responses  for  Case 
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4-24.  System  Matrices  for  Case  B-l,  Multipliers 

{1,1, 0,0, 1,1) 
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Fig.  4-27.  Control  Histories  for  Case  B-l,  Vertical  Translation  Maneuver  with  Multi- 


Time  Responses  for  Case  B-l,  Vertical  Translation  Maneuver  wi 


for  multipliers  results  in  the  matrices,  eigenvalues, 
and  time  histories  shown  in  Fig.  4-29,  Table  4-9,  and 
Figs.  4-30  and  4-31  for  g  =  1000.  An  interesting  result 
here  is  that  the  control  surface  deflections  (Fig.  4-30) 
both  "jumped"  to  a  maximum  value  and  then  quickly  returned 
to  equilibrium.  This  implies  some  type  of  new  trim  condi¬ 
tion  when  effecting  the  direct  lift  control  maneuver, 
which  is  unrealistic.  In  this  case  the  short  period 
approximation  is  physically  invalid  since  the  direct  lift 
control  maneuver  effects  a  forward  speed  change,  and  u 
is  assumed  to  be  zero  in  a  short  period  approximation. 

The  time  response  in  Fig.  4-31  shows  satisfactory  achieve¬ 
ment  of  the  final  commanded  values  but  with  a  great  deal 
of  output  interaction  in  the  first  two  seconds. 

In  an  attempt  to  reduce  this  interaction  the  eigen¬ 
vector  multipliers  are  changed,  with  the  best  case  found 
to  be  {1,1, 0,0, 1,1}.  This  is  the  same  as  the  best  case 
in  the  previous  test  case  and  produces  the  same  eigenvalue 
movements  as  shown  in  Fig.  4-24.  The  system  matrices, 
control  histories,  and  time  responses  for  g  =  1000  are 
shown  in  Figs.  4-32  through  4-34.  Fig.  4-34  (the  time 
responses)  still  indicates  interaction  between  the  outputs 
but  it  is  much  less  than  in  the  first  example. 


Case  B-3 

In  this  test  case  the  chosen  slow  eigenvalues  are 
{  (-4+j0) ;  (~4  +  j0) ;  (-2  +  j0)}  with  the  output  equations  given 
by  equation  (4-15) .  As  explained  in  Appendix  B,  generalized 
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4-29.  System  Matrices  for  Case  B-2,  Multipliers 

{1,0, 0,0, 1,0} 
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Table  4-9.  Eigenvalue  Movements  for  Case  B-2,  Multipliers 
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Fig.  4-32.  System  Matrices  for  Case  B-2,  Multipliers 

{1,1, 0,0, 1,1} 
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Fig.  4-33.  Control  Histories  for  Case 
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eigenvectors  are  not  generated  so  that  in  order  to  pro¬ 
duce  repeated  eigenvalues  (and  have  .  ,.Cn]  be 

invertible)  a  multiple  of  one  null  space  vector  is  used 
for  one  eigenvalue  a  multiple  of  the  other  vector  is  used 
to  produce  the  second  eigenvalue.  Testing  shows  that 
the  results  obtained  are  independent  of  the  multiple  of 
each  vector  used.  In  other  words,  the  same  results  are 
obtained  for  a  multiplier  choice  of  [a , 0 , 0 , b ,c ,d ]  where 
a  and  b  are  different  for  each  case  but  c  and  d  are  the 
same  (the  first  four  multipliers  correspond  to  the  re¬ 
peated  eigenvalues).  Choosing  (1 , 0 , 0 , 1 , 1 , 0 }  for  the  multi¬ 
pliers  yields  the  results  shown  in  Fig.  4-35,  Table  4-10, 
and  Figs.  4-36  through  4-43.  To  show  the  effect  of 
increasing  the  gain  parameter,  g,  Figs.  4-36  and  4-37 

show  the  control  and  output  histories,  respectively,  for 

T 

g  =  250  with  command  vector  v  =  II  0]  deg.  Figs.  4-38 

and  4-39  show  the  same  results  for  g  =  1000.  It  is  easily 

seen  that  the  response  of  y-^  =  8  becomes  tighter  with 

increasing  gain  while  the  amount  of  interaction  ( y ^ 

response)  changes  very  little.  The  same  results,  for 

T 

g  =  250  and  g  =  1000  but  with  v  =  [01]  deg,  are  shown  in 
Figs.  4-40  through  4-43.  Here  the  response  of  y2  =  Y 
remains  roughly  the  same  while  the  amount  of  interaction 
is  reduced  significantly.  These  results  show  that  the 
output  Y2  always  includes  a  significant  amount,  independent 
of  the  gain,  of  one  or  several  of  the  modes  defined  by  Z2 
(those  which  are  selected  as  transmission  zeros) . 
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4-35.  System  Matrices  for  Case  B-3,  Multipliers 

{1,0, 0,1, 1,0} 
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Table  4-10.  Eigenvalue  Movements  for  Case  B-3,  Multipliers 
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Control  Histories  for  Case  D-3,  Pitch  Pointing  Maneuver  with  Multipliers 
{1,0,0,1,1,01  (g  =  250) 
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Fig.  4-37.  Time  Responses  for  Case  B-3,  Pitch  Pointing  Maneuver  with  Multipliers 
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Case  B-4 


f 


This  case  is  identical  to  the  last  (same  eigenvalues 
and  eigenvectors)  except  that  the  output  equations  of 
equation  (4-16)  are  used  to  achieve  the  direct  lift  con¬ 
trol  mode.  The  system  matrices  are  shown  in  Fig.  4-44, 
and  the  eigenvalue  movements  are  the  same  as  in  Table 
4-10.  The  control  and  output  histories  for  g  =  1000 
are  shown  in  Figs.  4-45  and  4-46,  respectively.  The 
responses  clearly  show  a  significant  amount  of  interaction 
in  the  first  two  seconds.  To  improve  this,  the  eigen¬ 
vectors  would  have  to  be  modified,  but  this  was  not  done 
here. 


System  C  -  Full  Aircraft  Model  With  Control  Actuator 
Dynamics 

In  this  group  of  tests  the  full  longitudinal  air¬ 
craft  model  is  used  (that  is,  the  u  state  is  added) .  The 
state  equations  of  Fig.  3-5  were  implemented  originally, 
but  since  a  is  desired  in  the  output  equations,  these 
equations  are  transformed  to  replace  w  with  a  by  the 
relationship  of  equation  (3-2).  This  transformation 
results  in  the  following  partitioned  matrices. 
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4-44.  System  Matrices  for  Case  B-4,  Multipliers 
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It  is  important  to  note  here  that  if  the  control  actuator 
dynamics  had  not  been  modeled,  the  plant  would  not  have 
been  in  the  proper  form  for  implementation  (see  equation 
(2-1))  and  a  transformation  would  have  been  necessary. 
Therefore,  modeling  control  actuator  dynamics  guarantees 
that  the  plant  is  in  the  proper  form  for  use  of  this 
method . 

The  output  equations  required  for  pitch  pointing 
and  vertical  translation  maneuvers  are  given  by 

T  fo  0  0  1  [  0  0 

[C, JC,]x(t)  =  [6  y]  =  i  x(t)  (4-19) 

z  0-101100 


and  those  for  direct  lift  control  by 


[C  'C  ] =  [ y  a ] 1  =  1 

*  oiooloo 


0  -1  0  1  i  0  0 


(4-20) 


The  matrix  C ^  is  obviously  rank  deficient  in  both  cases, 
thus  again  requiring  extra  plant  measurements. 

For  these  tests,  four  slow  eigenvalues  had  to  be 
selected,  thus  requiring  eight  constant  multipliers  to 
define  the  closed-loop  eigenvectors. 
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Numerous  test  cases  were  completed  using  this  system. 
Since  all  the  results  are  similar,  only  one  test  case  is 
developed.  The  four  chosen  eigenvalues  are  {(-4±3j); 
(-4+jO) ;  (-3+jO)},  with  output  matrices  given  by  equation 

(4-19) .  The  first  selection  of  null  space  vector  multi¬ 
pliers  is  {1 , 0 , 0 , 0 , 1 , 0 , 1 , 0} .  The  system  matrices  are 
shown  in  Fig.  4-47,  and  the  movements  in  Table  4-11.  Note 
that  the  determinant  of  the  F2B2  matrix  is  very  small, 
therefore  causing  the  matrix  to  have  extremely  large 
terms  (relation  given  in  equation  (2-38)).  The  response 
for  this  system  is  totally  unrealizable.  Due  to  the  mag¬ 
nitude  of  the  response  no  plots  are  shown,  but  the  results 
for  g  =  1000  require  thousands  of  degrees  of  control  de¬ 
flections  and  have  output  response  overshoots  as  high  as 

T 

700  degrees  for  v  =  [10]  deg  and  1200  degrees  for  v  = 

T 

[0  1]  deg.  After  roughly  six  seconds  the  control  deflec¬ 
tions  do  settle  out  to  reasonable  values  and  the  outputs 
assume  their  commanded  values. 

For  multipliers  {1,10,0,0,0,1,0,1}  the  system 
matrices  are  shown  in  Fig.  4-48  and  eigenvalue  movements 
in  Table  4-12.  The  control  deflections  and  output 
responses  for  this  choice  are  still  unrealizable  but  are 
reduced  in  peak  magnitude  by  roughly  a  factor  of  ten.  No 
combination  of  multipliers  was  found  which  produced  a 
reasonable  response,  but  the  amount  of  testing  done  was 
far  from  exhaustive. 
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Table  4-11.  Eigenvalue  Movements  for  Case  C-l,  Multipliers  {1 
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Fig.  4-48.  System  Matrices  for  Case  C-l,  Multipliers 

{1,10,0,0,0,1,0,1} 
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Table  4-12.  Eigenvalue  Movements  for  Case  C-l,  Multipliers  {1,10, 


System  D  -  Robustness  Tests 


In  order  to  examine  how  well  this  method  behaves 
when  variations  in  the  system  plant  are  encountered,  a 
series  of  robustness  tests  are  performed.  Appendix  A-2 
shows  the  values  of  the  aircraft  stability  derivatives 
when  the  design  altitude  and  Mach  number  are  changed  from 
h  =  3000  ft,  M  =  0.60  to  h  =  10,000  ft,  M  =  0.70.  Also, 
in  the  test  cases  encorporating  control  actuator  dynamics, 
the  actuator  transfer  functions  are  also  changed  from 
that  of  equation  (3-3)  to 


&  _  15 

T~  s+15 


(4-21) 


To  demonstrate  robustness,  these  altered  plants  are  used 
with  the  controller  and  transducer  matrices  found  in  all 
of  the  previous  test  cases,  except  for  the  full  model 
test  cases  where  no  satisfactory  results  were  obtained. 

The  aircraft  plants  used  in  these  tests  are  found 
in  Appendix  A-2.  In  all  of  the  tests  for  robustness,  no 
noticeable  change  in  time  histories  occurs.  No  plots  are 
shown  because  they  appear  to  be  exact  duplicates  of  the 
design  condition  results. 


V.  Concl usions  and  Recommendations 

Nine  major  conclusions  are  drawn  from  the  results 
presented  in  this  study.  A  number  of  these  conclusions 
suggest  that  a  great  deal  more  study  is  necessary  before 
all  the  effects  of  the  various  aspects  of  this  method  are 
fully  understood.  The  conclusions  are  presented  in  a 
designer's  order,  that  is,  they  try  to  answer  the  questions 
that  would  arise  as  the  controller  design  process  proceeds. 

The  major  effect  of  the  selection  of  the  slow  eigen¬ 
values,  beside  the  fact  that  they  insure  stability  at 
high  gain  when  chosen  to  be  in  the  left-half  plane,  is 
that  their  selection  changes  the  magnitude  of  the  F2B2 
matrix.  For  any  given  set  of  selected  slow  eigenvalues, 
the  determinant  of  ^oes  not  change  for  different 

eigenvector  choices.  Since  Kq  is  found  by  inverting  F2B2' 
the  larger  the  value  of  the  determinant,  the  smaller  the 

terms  in  K  and  K.  will  be.  The  size  of  the  terms  in  K 
o  1  o 

has  a  major  effect  on  the  response,  as  seen  in  many  of  the 
test  cases. 

The  selection  of  the  slow  eigenvectors  directly 
affects  the  transient  response  of  the  outputs.  The  effect 
of  the  choice  of  the  eigenvectors  is  best  seen  in  Cases 
A-l,  A-3,  and  B-l.  Different  eigenvectors  produce  dif¬ 
ferent  M  matrices  and  therefore  different  measurement 
matrices,  and  F2.  Although  the  determinant  of  F2B2 
remains  constant  for  any  choice  of  eigenvectors,  the 


elements  of  ^2^2  c^an9e  anc^  thus  the  elements  of  Kq  and 
change.  These  changes  in  the  controller  matrices  radically 
alter  the  transient  response.  It  has  been  observed  in  all 
the  test  cases  which  showed  different  eigenvector  choices 
that  a  system  with  larger  values  of  Kq  and  has  a  poorer 
response  than  one  with  smaller  values  in  these  matrices. 

The  choice  of  the  elements  in  the  E  matrix  was 
briefly  examined  in  Case  A-6.  The  only  effect  of  changing 
the  values  of  E  was  to  raise  or  lower  the  value  of  g  by 
the  inverse  of  the  amount  by  which  o  was  changed.  That 
is,  if  E  is  decreased  by  a  factor  of  ten,  the  responses 
are  identical  except  that  the  value  of  g  must  be  increased 
by  a  factor  of  ten.  None  of  the  test  cases  examined  the 
effect  of  choosing  f  a ^  or  changing  one  value  of  0 
while  holding  the  other  constant. 

Changing  the  ratio  between  Kq  and  has  a  small 
but  noticeable  effect  on  the  response.  The  ratio  of 
to  Kq  defines  the  first  set  of  slow  eigenvalues  (see 
equation  (2-39)).  Test  cases  A-l  and  A- 3  along  with  A- 2 
and  A-4  are  identical  except  that  in  A-l  and  A-2,  Kq  = 
and  in  A-3  and  A-4,  =  2Kq.  The  effect  of  this  change 

is  to  slightly  decrease  the  final  value  overshoots  and 
settling  times,  which  suggests  the  first  set  of  slow  eigen¬ 
values  are  slightly  visible  in  the  response. 

The  gain  parameter,  g,  has  the  single  most  pro¬ 
nounced  effect  on  the  controller  performance.  In  all 
the  aircraft  plants  used  in  this  study,  the  open-loop 
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response  is  unstable  due  to  a  right-half  plane  pole.  As 
g  increases,  one  of  the  integrator  poles  (which  is  at  the 
origin  in  the  open-loop)  moves  into  the  right-half  plane, 
meets  with  the  unstable  pole,  and  they  then  split  off  the 
real  axis  and  move  back  toward  the  left-half  plane.  At 
some  value  of  g  they  enter  the  left-half  plane  and  the 
system  becomes  stable.  As  g  is  increased  further,  the 
system  poles  approach  their  asymptotic  values.  However, 
as  g  is  made  very  large,  the  terms  in  the  fast  state 
differential  equations  become  very  large  and  integrating 
them  becomes  both  costly  and  inaccurate.  In  several  of 
the  cases  in  this  report,  very  large  gains  were  tried, 
which  resulted  in  the  computed  response  degrading  from 
the  results  shown.  This  is  possibly  the  result  of  trying 
to  numerically  integrate  a  system  of  equations  that  in¬ 
cludes  "stiff"  modes. 

In  all  of  the  B  test  cases  and  in  the  full  model 
tests,  control  actuator  dynamics  are  modeled  in  the  air¬ 
craft  state  equations.  The  responses  in  the  A  test  cases 
are  too  fast — modeling  actuator  dynamics  in  the  B  cases 
slowed  the  responses  to  an  acceptable  level.  In  the  full 
model  cases  the  basic  plant  (without  actuator  dynamics) 
is  not  in  the  form  required  to  apply  Porter's  method.  By 
introducing  actuator  dynamics,  one  can  avoid  the  necessity 
of  transforming  the  state  equations  since  the  plant  is 
automatically  put  into  the  proper  form. 


The  robustness  tests  demonstrate  the  most  attrac¬ 


tive  feature  of  this  method.  In  all  cases,  when  the  flight 
condition  is  slightly  changed  and  the  model  of  the  control 
actuators  altered,  the  original  controller  and  transducer 
matrices  still  produced  the  desired  flight  modes.  Not 
only  are  these  modes  satisfactorily  effected,  the  response 
is  practically  identical  to  that  of  the  original  system. 

Unfortunately,  not  all  desired  outputs  can  be  tracked 
with  this  method.  For  example,  when  pitch  rate  q,  or  nor¬ 
mal  acceleration  A^,  is  a  desired  output,  it  is  not 
possible  to  find  an  M  matrix  which  satisfies  the  necessary 
conditions  (that  the  second  set  of  slow  eigenvalues  are 
all  in  the  left-half  plane  and  that  the  matrix  F2  has 
full  rank) . 

The  most  accurate  outputs  to  describe  the  three 
modes  in  this  study  are  not  6,  y,  and  a,  but  q,  A^,  and  a. 
Since  q  and  AN  outputs  yield  uncontrollable  systems,  they 
were  replaced  with  6  and  y  as  approximations.  The  reason 
for  this  failure  is  that  the  system  matrix  for  these  out¬ 
puts  is  not  fully  controllable,  so  the  method  will  not 
work.  These  cases  (when  q  or  A^  was  used)  represented  the 
only  ones  in  which  F£  was  not  full  rank. 

The  final  conclusion  deals  with  the  full  model  of 
the  aircraft.  For  some  reason,  it  is  extremely  difficult 
in  this  case  to  choose  the  eigenvalues,  the  eigenvectors, 
the  £  matrix,  and  a  ratio  of  and  that  produces  a 
satisfactory  response.  The  only  major  trend  noticeable 
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here  is  that  the  response  is  drastically  worse  when  the 
controller  matrices  contain  very  large  terms. 


Recommendations 

The  full  effects  of  all  the  design  parameters  in 
this  method  still  require  a  great  deal  of  investigation. 
Particularly,  the  effect  of  the  chosen  eigenvalues  and 
eigenvectors  on  the  response  still  requires  more  study 
since  the  results  clearly  indicate  they  do  appear  in  the 
closed-loop  system.  Also,  the  effects  of  different 
choices  of  I  should  be  examined  further.  The  problems 
with  the  full  model  cases  require  particular  attention. 

It  would  be  highly  beneficial  to  the  advancement  of  this 
method  if  a  standard  way  to  choose  the  various  design 
parameters  to  minimize  the  amount  of  commanded  value  over¬ 
shoot  and  interaction  could  be  devised.  Since  choosing 
different  eigenvectors  greatly  affects  the  transient 
responses,  it  is  possible  that  the  work  being  done  on 
entire  eigenstructure  assignment  problems  may  be  appli¬ 
cable  to  this  method. 
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Appendix  A-l.  Data  for  the  AFTI-16  Aircraft 


Mass  =  653  slugs 

C  =  11.32  ft 

S  =  300  ft2 

Iy  =  53876  slug/ft2 

b  =  30  ft 

Flight  Condition  Data 

M  =  0.60 

a  =  1104  ft/sec2 

h  =  3,000  ft 

U  =  663  ft/sec 

o  =  0.002177 

slug/ft2 

q  =  478  lb/ft2 
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Appendix  A-2.  Robustness  Test  Models 


Physical  Description  of  the  Aircraft 


Mass  =  653  slugs 

C  =  11.32  ft 

S  =  300  ft2 

Iy  =  53876  slug/ft 

b  =  30  ft 

Flight  Condition  Data 

M  =  0.70 

a  =  1077  f t-sec~2 

h  =  10,000  ft 

U  =  754  ft-sec  2 

p  =  0.001756  slug-ft-3 

q  =  499  lb-ft-2 

Non-Dimensional  Stability 

Derivatives 

Assumed  to  remain  the  same  as  in  Appendi 

to  small  flight  condition 
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Appendix  B.  Discussion  of  the  Usage  of 
Null  Space  Vectors  from  a  Singular  Value 
Decomposition  to  Form  Closed-Loop  Eigenvectors 


Real,  Ncn-repeated  Eigenvalues 

When  all  the  chosen  slow  eigenvalues  are  real,  and 
none  are  repeated,  the  null  space  o  f  [A.I-A^,  Ax  2 3  cor¬ 
responding  to  each  assigned  eigenvalue  contains  £  vectors, 
each  of  which  has  length  n+£  (£  is  the  number  of  inputs, 
n  is  the  number  of  states).  Therefore,  each  set  of  £ 
vectors  must  be  linearly  combined  (requiring  £  constant 
multipliers)  to  form  each  eigenvector.  The  only  restric¬ 
tion  placed  upon  the  multipliers  in  this  case  is  that  the 
eigenvectors  produced  are  linearly  independent. 


Complex  Eigenvalues 

The  selection  of  complex  eigenvalues  results  in 
complex  eigenvectors.  Since  the  singular  value  decomposi¬ 
tion  computer  routine  works  only  with  real  numbers,  the 
complex  eigenvectors  must  be  selected  from  vectors  lying 
in  the  null  space  of  (Ref  6) : 


0 


(B-l ) 


where  the  subscripts  r  and  i  refer  to  real  and  imaginary, 
respectively.  Using  these  vectors,  the  feedback  gain 
matrix  S  is  found  from 
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f  •  /  W  -  / 

X1  ^ 


’— n  ^ 


(B-2 ) 


In  this  case,  however,  21  vectors  lie  in  the  null  space 
of  the  matrix  in  (B-l)  for  each  complex  eigenvalue.  This 
is  due  to  the  fact  that  complex  eigenvalues  come  in  pairs, 
with  half  of  the  vectors  corresponding  to  the  eigenvalue 
with  the  positive  imaginary  part  and  the  other  half  to  the 
eigenvalue  with  the  negative  imaginary  part.  Furthermore, 
each  of  these  21  vectors  has  2(n+£)  elements,  corresponding 
to  the  real  and  imaginary  part  of  each.  Therefore,  the 
vectors  should  be  combined  with  21  multipliers  and  then 
are  separated  to  conform  with  equation  (B-2) .  In  this 
study,  however,  half  of  the  null  space  vectors  are  given 
zero  multipliers  and  the  £  remaining  vectors  are  combined 
and  separated  to  form  the  closed-loop  eigenvectors.  This 
is  found  to  correctly  reproduce  the  chosen  eigenvalues 
when  the  corresponding  closed- loop  A  matrix  is  formed,  so 
its  use  is  justified. 

Repeated  Eigenvalues 

If  repeated  eigenvalues  are  desired,  two  approaches 
are  possible  depending  on  how  many  repeated  values  are 
desired.  If  the  number  of  desired  repeated  eigenvalues 
is  less  than  or  equal  to  the  dimension  of  the  null  space 
(£) ,  then  it  is  possible  to  choose  linearly  independent 
eigenvectors  from  these  spaces.  In  this  case,  no  exten¬ 
sion  to  the  previously  developed  theory  is  necessary, 
but  the  designer  is  more  restricted  in  eigenvector 
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selection.  One  way  to  avoid  choosing  linearly  dependent 
eigenvectors  (the  one  used  in  this  study)  is  to  choose 
each  eigenvector  as  a  multiple  of  one  of  the  null  space 
vectors . 

To  assign  more  than  l  repeated  eigenvalues  or  to 
increase  the  flexibility  in  eigenvector  selection  for  £ 
or  less  repetitions,  a  string  of  generalized  eigenvectors 
must  be  generated.  Since  it  is  not  used  in  this  study, 
the  reader  is  referred  to  Ref  6  for  further  details. 
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Appendix  C.  Main  Computer  Program  and  Example 
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